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The recent discovery of superconductivity at 190 K in highly compressed H 2 S is spectacular not 
only because it sets a record high critical temperature, but because it does so in a material that 
appears to be, and we argue here that it is, a conventional strong-coupling BCS superconductor. 
Intriguingly, superconductivity in the observed pressure and temperature range was predicted the¬ 
oretically in a similar compound H 3 S. Several important questions about this remarkable result, 
however, are left unanswered: (1) Does the stoichiometry of the superconducting compound differ 
from the nominal composition, and could it be the predicted H 3 S compound? (2) Is the physical 
origin of the anomalously high critical temperature related only to the high H phonon frequencies, 
or does strong electron-ion coupling play a role? We show that at experimentally relevant pressures 
H 2 S is unstable, decomposing into H 3 S and S, and that H 3 S has a record high due to its cova¬ 
lent bonds driven metallic. The main reason for this extraordinarily high in H 3 S as compared 
with MgB 2 , another compound with a similar superconductivity mechanism, is the high vibrational 
frequency of the much lighter H atoms. 


Recently reported superconductivity at 190 K in com¬ 
pressed H 2 S [U has been arguably the biggest discovery in 
the field since the superconducting cuprates nearly thirty 
years earlier. Superconductivity in a related compound, 
H 3 S, in a similar range of pressure with very nearly the 
same critical temperature was predicted theoretically [5] 
about a year earlier. In that theoretical paper, direct ab 
initio calculations yielded high phonon frequencies, giv¬ 
ing the logarithmic average (the pre-factor in the equa¬ 
tion for Tc) on the order of 1100-1300 K and a coupling 
constant A larger than 2, combining to give between 
191-204 K. However, a microscopic understanding of why 
this particular material features such a strong coupling is 
still missing, as is an explanation of the discrepancy be¬ 
tween experimental and theoretical stoichiometries. As 
we show here, the answer lies in the stability of the H„S 
series of compounds. At high pressure, the phase di¬ 
agram favors decomposition of H 2 S into H 3 S and pure 
S. The mechanism of superconductivity can be traced 
to the strongly covalent metallic nature of H 3 S along 
with high phonon frequences, similar to another con¬ 
ventional (i.e. phonon-driven) supercondnctor with what 
now seems only a relatively high T^, MgB 2 [3]. 

The discovery in 2001 of phonon-driven superconduc¬ 
tivity at 39 K in MgB 2 not only set a record high Tc 
for a conventional phonon-mediated mechanism, which 
just 30 years back was widely believed to be limited 
to < 25 K, but also introduced a completely new con¬ 
cept in the theory of superconducting materials, dnbbed 
“doped covalent bonds” by W. E. Pickett and his col¬ 
laborators [1]. The essence of this concept is that bond¬ 
ing and antibonding states in a covalent system are very 
sensitive to hopping integrals and thus to ionic positions, 
which makes them strongly coupled to the corresponding 
phonons. However, in essentially all covalent systems, the 
corresponding states are removed from the Fermi surface 
and thus irrelevant to superconductivity. Moreover, even 


when it is possible to dope a covalent insulator, such as 
diamond, it costs a tremendous amount of energy and 
therefore results in only very small doping levels. MgB 2 
is different in two ways: First, it sports metallic tt bands 
in addition to the strongly covalent a bands; second, the 
a bands are 2D, and therefore even a small carrier concen¬ 
tration in these bands creates a sizeable density of states 
(DOS), that is, a substantial covalent metallicity. On the 
contrary, in diamond all the bands are 3D, In contrast, all 
the bands are 3D in diamond, and small doping causes a 
similarly small DOS, and thus a low critical temperature 
T 

c- 

While the general concept that hard phonons in H- 
rich materials might make them good superconductors is 
not new, so far high-Tc superconductivity has been elu¬ 
sive, primarily because such hard phonons normally do 
not produce large coupling constants. As V. L. Ginzburg 
wrote in 1977, “... in many already-known materials, the 
Debye temperature is very large, ~ 10^ K, and low is 
related to small coupling constants ... In view of this, 
attention is attracted to various hydrogen-rich materials 
under high pressure.” [5] The strong electron-phonon cou¬ 
pling inherent to a covalent metal such as MgB 2 suggests 
a way to circumvent this problem, by combining met¬ 
allized covalent bonds with lighter elements and higher 
phonon frequencies. We argue that H 3 S applies precisely 
this recipe to achieve its record critical temperature and 
that the physics of superconductivity, and, to some ex¬ 
tent, even numerics, in H 3 S is extremely similar to MgB 2 
with the only qualitative difference being the factor of 11 
between the masses of H and B. In plain English, H 3 S is 
like MgB 2 , but lighter. 

To gauge the level of electron-phonon coupling, we 
note that the same electron-phonon interactions that 
contribute to phonon-mediated superconductivity also 
manifest themselves as the screening of bond-stretching 
force constants and softening of vibrations; this soften- 
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ing can therefore be taken as a proxy for the strength 
of the coupling. Recalling the case of MgB 2 , we observe 
that it has covalent-bond stretching phonons at about 
500 cm“^ near the zone center, while in its sister com¬ 
pound AIB 2 , which has its cr band Fermi level and un¬ 
coupled from conducting electrons, those phonons are as 
hard as 900 cm“^ [S]. A softening of similar magnitude 
occurs in solid H 3 S as compared with vibrations in S-H 
containing thiol molecules. The latter have frequencies 
of about 2500 cm“^, while the calculated phonon fre¬ 
quencies in H 3 S are roughly ~ 1600 cm“^ [5]. This is a 
very large softening indicating a very large coupling, even 
though it may not initially perceived as such because the 
bare frequency involving the light H atoms is so high. 

In any theoretical analysis of a new material, it is im¬ 
portant to establish the stoichiometry and the crystal 
structure of the compound of interest. For instance, a 
recently suggested theory [7] is based upon the H 2 S com¬ 
position, which, as we show below, is nearly certainly 
not the composition that supports superconductivity in 
the experiment. While the experiments showing super¬ 
conductivity at 190 K started by compressing H 2 S, they 
also showed the formation of pure S, suggesting that the 
material that actually superconducts is likely to be H 
enriched [1]. The composition that Duan et al. [2] stud¬ 
ied theoretically, H 3 S, is consistent with this observation, 
but they did not consider the full range of compositions 
in the phase diagram. Previous theoretical publications 
have verified that the proposed high pressure phases for 
H 3 S in Ref. [2] and H 2 S in Ref. [ 8 ] are stable against decom¬ 
position into H and S. However, lack of decomposition 
into elemental species is not a particularly stringent test, 
and stability against separation into other phases, e.g. 
H 2 S into S and H 3 S, which is important for understand¬ 
ing the relevance of any calculations to the experiments 
in Ref. [TJ has not been investigated until now. 

We begin by checking the stability of H„S compounds 
with respect to decomposition into other phases by cal¬ 
culating their zero temperature enthalpy H — E + PV as 
a function of composition using density functional theory 
calculations. We used the VASP density functional the¬ 
ory software [9] with the Perdew-Burke-Enzerhoff gen¬ 
eralized gradient approximation m. and a projector- 
augmented waves basis [nnn] with a 1000 eV plane wave 
cutoff. The calculations used 16 x 16 x 16 Monkhorst- 
Pack fc-points for cells containing a single formula unit; 
correspondingly reduced A:-point grids were used for the 
larger cells. Geometries were relaxed under an applied 
pressure using the conjugate-gradient algorithm applied 
to both unit cell size and shape, as well as atomic posi¬ 
tions, until the residual forces were less than 0.01 eV/A. 
Since the stability of each composition depends on the 
enthalpy of the lowest enthalpy structure at that com¬ 
position, we considered many structures at compositions 
ranging from pure S, to H„S for n = 1 — 6 , to pure 
H. The initial structures for our relaxation procedure 


came from previously published experimental and com¬ 
putational studies [aiaiiaHii], manual modifications of 
these published structures, and from a simple version 
of the random structure search method [H]. The final 
relaxed structures are listed in the Supplemental Mate¬ 
rial [17] . 

The zero temperature formation enthalpy of HjjSi-a,, 
AH{x) = H{x) — xHyi — (1 — x)Hs as a function of x, is 
plotted in Fig. we see that at P = 200 GPa, H 3 S 
in the previously proposed body-centered-cubic (bcc) 
/mSm structure is stable with respect to decomposition 
into any of the other calculated structures. At 200 GPa 
the i?3m structure has cell and internal parameters that 
make it identical to the /m3m, but below about 150 GPa 
it begins to continuously evolve to a distinct, lower sym¬ 
metry structure m- All simulated H 2 S structures, on 
the other hand, are unstable with respect to decomposi¬ 
tion into S and H 3 S by at least 120 meV/formula unit. 
All structures with higher H concentration lie above the 
convex hull, and are therefore also unstable with respect 
to decomposition into H 3 S and H. In fact, all such struc¬ 
tures include H atom pairs with distances similar to that 
of the H 2 molecule showing the tendency to phase 
separation and release of pure H at these compositions. 
Note that small errors in the enthalpies of the end points 
(pure S and pure H), if for example the true structures are 
a bit lower in enthalpy than the ones we considered, will 
not change these conclusions with respect to the stability 
of H 3 S and instability of all other compositions. How¬ 
ever, it is in principle possible that a sufficiently lower 
enthalpy H structure could make all intermediate com¬ 
positions unstable, or that a H 2 S structure with enthalpy 
outside of the convex hull, which has not been consid¬ 
ered here, could exist. We think that the existence of 
such a lower enthalpy structure is very unlikely for pure 
H, which has been studied extensively [15]. For H 2 S, 
which is chemically reasonable and has not been previ¬ 
ously studied extensively under such high pressure, we 
think that the failure of our thorough intuition-guided 
and random searches to find such a lower enthalpy struc¬ 
ture makes its existence highly unlikely. Finally, higher H 
concentration phases show a tendency for formation of H 2 
molecules consistent with their predicted instability with 
respect to decomposition into pure H and H 3 S. Reducing 
the pressure to 150 GPa does not change the enthalpies of 
the low-lying structures significantly. We therefore con¬ 
clude that under experimentally relevant pressures the 
starting material H 2 S decomposes as 

3 H 2 S 2 H 3 S + S, (1) 

with no other products. 

An important issue we have not addressed so far is 
the reliability of projector-augmented wave (PAW) cal¬ 
culations at such compressions since, generally speaking, 
the available PAWs were designed and tested for much 
larger interatomic separations. To this end, we compared 
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FIG. 1. (color online) Formation enthalpy of as a 

function of H concentration at P = 150 GPa (red x’s) and 
200 GPa (black squares). Upper axis tick marks indicate com¬ 
positions equivalent to HnS for integer n. All compounds 
above the convex hull (lines) are unstable with respect to de¬ 
composition into the two adjacent phases on the convex hull. 
At both pressures the only stable compound is H 3 S. Pure S 
structures are taken from Refs. and m and pure H from 
Ref. [151 



FIG. 2. (color online) Band structure of the H 3 S 7m3m struc¬ 
ture at P = 200 GPa, calculated using FPLO. Weights of the 
three most important atomic orbitals are shown with the sym¬ 
bols of the corresponding sizes. The large contribution of H 
orbitals to the S-s derived states at the bottom of the band 
indicates strong covalency. 


the energy difference of the thermodynamic reaction Q 
with the energy difference calculated with an all-electron 
method, full-potential local orbitals (FPLO) [TS]. We 
found a formation energy of 105 meV/formula unit, 
within 10% of the VASP value, confirming that the latter 
calculations are reliable. 

Having established that the material exhibiting super¬ 
conductivity at 190 K at P ~ 200 GPa is most likely 
the Im3m bcc-like structure of H 3 S identified in Ref. HI 
we analyze its electronic structure and superconductivity. 
We do not aim at recalculating the exact numbers for the 
coupling constants and phonon frequencies, but rather at 
gaining insight into why the calculations of Ref. [2] pro¬ 
duced such a large coupling and T^. 

In Fig. we show our calculated band structure, 
which agrees with Duan et al. The Fermi surface for 
the one band (out of the five that cross the Fermi 
level) that contributes the overwhelming majority of 
the density of states, colored by the Fermi velocity, is 
shown in Fig. The calculated average Fermi veloc¬ 
ity, defined as vp = -\/Sk<J(Pk — Ep)v'^/Ti\^5{E\^ — Ep) 
is 0.25 X 10® cm/sec. This allows us to address an¬ 
other question: given the large Tc, would the coher¬ 
ence length be long enough for the standard Eliash- 
berg theory to be applicable? We know that in high- 
Tc cuprates this is nearly the case, which was argued to 
have important theoretical implications in terms of Bose- 
condensation of local pairs rather than BCS long-range 
coherence |19j. For H 3 S, we can estimate the zero temper¬ 
ature gap parameter, A(0), using Carbotte’s formula PO] . 
A(0) = 1.76fcBTe[l + 12.5(T,/wiog)2log(wiog/2r,)]. Us¬ 
ing the numbers from Ref. HI we get A(0) « 40 meV. 



FIG. 3. The main pocket of the Fermi surface of H 3 S at 
P = 200 GPa, colored according to the local Fermi velocity 
(the scale is in arbitrary units). Note heavy bands near the 
van Hove singularities. Two small hole pockets near F and 
an electron pocket near H contribute little to the total DOS, 
and are omitted from the picture. 


Using the standard expression for the clean limit coher¬ 
ence length, ^ = hvp/TrA{0), we find f ~ 40 A, much 
larger than the interatomic distance. 

Analyzing the characters of the wave functions as in 
FigH we observe that the bands at the Fermi level are 
formed nearly exclusively by seven orbitals: sulfur 3s, 
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sulfur 3px,y,z, and the three Is orbitals of the hydrogens, 
each displaced along x, y, or z from its nearest neighbor 
sulfur. The S d—states are located more than 15 eV above 
the Fermi and can be safely neglected. In the following 
we will denote the three hydrogen Is orbitals hy i = 
X, y, z, and the three sulfur p orbitals hy I = X, Y, Z. 
The nearest neighbor S-H Hamiltonian is 

H,, = taC,, (2) 

Hit = hSi, (3) 

where ta,b are the sstr and psa S-H hoppings, respec¬ 

tively, Ci = 2cos{kia/2), and Si = 2sm{kia/2). The 
H-H nearest neighbor hopping, despite this distance be¬ 
ing the same as the S-H one, is much smaller, due to the 
difference in radii of the H Is and S 3s orbitals. The H-H 
Hamiltonian can be written down as 

Hxy = tcexp[i{kx - ky)a/2]Cz, (4) 

etc. Using FPLO to construct these Wannier functions 
and their corresponding tight-binding Hamiltonian, we 
find ta = —4.2 eV, 4 = —5.2 eV, and tc = —2.7 eV. 
The onsite energies are: = — 8.6 eV, E^p = —1.3 eV, 

and F1 is(h) = —5.0 eV {Ep set to zero). We see that the 
hydrogen s levels lie between the two sulfur levels, and 
create strong covalent bonds with both. The calculated 
bonding-antibonding splitting between sulfur p and hy¬ 
drogen s states at the P point (fc = { 7 r/ 4 o, 7 r/ 4 a, 7 r/ 4 a}) 
is 2A « 25 eV. However, the same splitting, by sym¬ 
metry, is zero at the F point, so, despite very strong 
covalency, this bond remains metallic, with the DOS 
iV(0) ^ 0.6 states/eV. 

It is now instructive to compare these parameters 
with those in MgB 2 , keeping only bond-stretching boron 
phonons in MgB 2 and bond-stretching sulfur phonons in 
H 3 S; in both cases these contributed about 70% to the 
total coupling (it is worth noting that a 30% contribution 
to the coupling constants does not imply a comparable 
contribution to Ty, in fact, since wiog = exclud¬ 

ing sulfur phonons would lead to very small changes in 
Tc). In order to do so, we use the well known qualitative 
relations between the electron-phonon coupling constant 
A, its electronic part (also called Hopfield factor) ry, and 
the average force constant $. For the purpose of the en¬ 
suing discussion it is enough to know that the Hopfield 
factor characterizes the electron-ion interaction and de¬ 
pends only on electronic properties, such as the DOS and 
the ionic potential, but not on phonon frequencies, while 
$ is defined in such a way m that it represents a com¬ 
bination of the derivatives of total energy with respect to 
ionic coordinates (force constants) and thus carries infor¬ 
mation about the phonon spectrum, but does not directly 
depend on the electronic energies and wave functions. On 
a semiquantitative level [52], 


A « 77/4> (5) 

$ = $0 - 277. (6) 

Here $0 represents unscreened force constants that do 
not account for electron-phonon coupling. Note that if 
there is one dominant phonon mode, $o/4> « 
where uj and ujq are the screened and unscreened frequen¬ 
cies of this mode, respectively. Using the bond-stretching 
mode frequencies for AIB 2 and MgB 2 to represent un¬ 
screened and screened bonds, respectively, as mentioned 
above, we get 

$o/4> « (900/500)2 _ 3 (- 7 ^ 


which qualitatively agrees with the accepted values for 
MgB 2 [3], Aero- =0.78-1.02. This gives 

VMgB.i « 1 X 11 X (500 cm“^)2 (9) 

~ 2.75 X 10®TO//-cm-^. (10) 

If we assume the same Hopfield factor 77 for H 3 S, we get 
A « 2.75 X 10®/13002 « 1.6 (estimating the average fre¬ 
quency of the bond-stretching modes in H 3 S from Fig. 5 
in Ref. [2] as 40 THz), which is in the right ballpark com¬ 
pared to the value A = 2.19 reported in Ref. |2| Note 
that the logarithmic frequency in that paper was cal¬ 
culated to be 930 cm“^, suggesting that phonons softer 
than 40 THz contribute to total coupling. Had we taken 
UJ = IlOO cm“^ = 35 THz instead of 40 THz, we would 
have obtained the value calculated in Ref. |2| 

As a consistency check, let us now see whether 
this estimate implies a plausible number for $ 0 - Us¬ 
ing w = 1300 cm“^ and A = 1.6, we find ujq « 
1300-\/l + 2 X 1.6 = 2660 cm“^, consistent with the vi- 
bron frequencies in S-H molecules, which have a large 
gap and are unscreened by definition. The internal con¬ 
sistency of our analysis and its consistency with experi¬ 
ment confirm our main conclusions below. 

In conclusion, at a pressure of P = 200 GPa, H 2 S gains 
nearly 40 meV per atom by decomposing into elemental 
sulfur and H 3 S in the Imim structure. The physical 
mechanism underlying the high-temperature supercon¬ 
ductivity of H 3 S at that pressure is very similar to that 
in MgB 2 : metallization of covalent bonds. The main dif¬ 
ference from MgB 2 is that the hydrogen mass is 11 times 
smaller than the mass of boron, resulting in a 3.5 times 
larger prefactor. 

We acknowledge useful discussions with D. A. Papa- 
constantopoulos and M. Calandra. 
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